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Abstract: 

Aerosol size distributions were measured by a Scanning Mobility Particle Sizer (SMPS) 
onboard the CIRPAS Twin Otter aircraft during 16 flights at the Southern Great Plains 
(SGP) site in northern central Oklahoma as part of the Aerosol Intensive Operation 
period in May, 2003. During the same period a second SMPS was deployed at a surface 
station and provided continuous measurements. Combined with trace gas measurements 
at the SGP site and back- trajectory analysis, the aerosol size distributions provided 
insights into the sources of aerosols observed at the SGP site. High particle 
concentrations, observed mostly during daytime, were well correlated with the sulfur 
dioxide (SO 2 ) mixing ratios, suggesting nucleation involving sulfuric acid is likely the 
main source of newly formed particles at the SGP. Aerosols within plumes originating 
from wildfires in Central America were measured at the surface site. Vertically compact 
aerosol layers, which can be traced back to forest fires in East Asia, were intercepted at 
altitudes over 3000 meters. Analyses of size dependent particle growth rates for four 
periods during which high cloud coverage was observed indicate growth dominated by 
volume controlled reactions. Sulfate accounts for 50% to 72% of the increase in aerosol 
volume concentration; the rest of the volume concentration increase was likely due to 
secondary organic species. The growth law analyses and meteorological conditions 
indicate that the sulfate was produced mainly through aqueous oxidation of SO 2 in clouds 
droplets and hydrated aerosol particles. 
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1. Introduction 


2 Atmospheric aerosols affect the global energy budget by scattering and absorbing 

3 sunlight (direct effect) and by changing the microphysical structure, lifetime, and 

4 coverage of clouds (indirect effect). While it is widely accepted that aerosols could have 

5 significant impact on global climate, at present the magnitudes of these effects are poorly 

6 understood. Unlike green-house gases, whose radiative forcing can be calculated with 

7 high accuracy, there are substantial uncertainties associated with predicted aerosol 

8 radiative forcing, especially the indirect forcing resulting from the interaction of aerosols 

9 and clouds. The Intergovernmental Panel on Climate Change (IPCC) estimated that the 

10 direct and indirect effects of aerosols remain the most uncertain components in forcing of 

11 climate change over the industrial period [IPCC, 2001], The large uncertainties in 

12 aerosol direct and indirect effects are combinations of the present fragmented 

13 understanding of the interactions among aerosol, radiation, and clouds, and the large 

14 uncertainties in knowledge of aerosol properties and global distributions. Compared to 

15 green-house gases such as CO 2 , the lifetime of atmospheric aerosol is very short, about a 

16 week. As a result, the temporal and spatial distributions of atmospheric aerosol and its 

17 properties are inhomogeneous and highly variable. Predictions of future aerosol effects 

18 require detailed understanding of aerosol properties and geographical distributions, as 

19 well as of the processes that control the evolution of aerosols. Among aerosol properties, 

20 aerosol size distribution is one of the most important parameters needed for assessing 

21 aerosol climate effects. Aerosol size distribution influences climate forcing through its 

22 effects on light scattering and modification of cloud microphysical properties affecting 

23 albedo and persistence. 
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1 The Southern Great Plains (SGP) Climate Research Facility is the first and largest 

2 of such facilities established by the Atmospheric Radiation Measurement (ARM) 

3 Program. The SGP Climate Research Facility covers an area of -142,000 km 2 in North 

4 Central Oklahoma and South Central Kansas. The central facility at SGP is located at 

5 36.60°N and 97.48°W, with an altitude of 319 m. Since 1992, continuous measurements 

6 have been carried out, which provide a long time record of aerosol, cloud, radiation, and 

7 atmospheric properties at the SGP. These measurements have been used in many aerosol 

8 studies including quantification of the aerosol indirect effect [ Feingold et al., 2003 and 

9 Kim et al., 2003], studies of vertical aerosol extinction and backscattering [ Turner et al., 

10 2001, 2002], and the application of surface measurements to retrieve aerosol optical 

11 depth [ Bergin et al., 2000]. From May 5 to 31, 2003, an Aerosol Intensive Operation 

12 Period (IOP) was conducted at the SGP [Ferrare et al., 2005a]. During the IOP, 

13 additional instalments were deployed both at surface sites and onboard a research 

14 aircraft. The additional measurements included sub-micrometer aerosol size distribution, 

15 hygroscopicity, chemical composition, and optical properties. In addition to aerosol 

16 measurements, the mixing ratios of relevant trace gases were also monitored at a NASA 

17 mobility facility during the IOP. The simultaneous measurements of the wide range of 

18 aerosol properties not only provided detailed characterizations of aerosols during spring 

19 time, but also allowed studies of aerosol evolution processes at the SGP. 

20 In this paper, we present both continuous records at surface and vertical profiles 

21 of aerosol size distribution measured during the IOP at the SGP site. Based on back- 

22 trajectory analyses, the air masses that arrived at the SGP site were grouped into 4 

23 classes. For each class, the source and the evolution process of the aerosol are discussed 
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using simultaneous aerosol and trace gas measurements. During the IOP, plumes from 

2 distant fire activities were detected both at the surface and in elevated layers at altitudes 

3 over 3000 m. Aerosols within these plumes showed high scattering coefficients and large 

4 accumulation-mode diameters and concentrations. The characteristics of the aerosols 

5 observed within these plumes are presented. We also carried out detailed particle growth 

6 law analyses for four cases using the continuous measurements of aerosol size 

7 distribution at the surface. The results of the growth law analyses provide insights into 

8 the aerosol evolution processes at the SGP. 

9 

10 2. Measurements. 

1 1 2.1 Aerosol size distribution measurements. 

12 Among the instruments deployed during the IOP were two Scanning Mobility 

13 Particle Sizers (SMPS). One of the two SMPS was operated onboard a Twin Otter 

14 aircraft to characterize the vertical profiles and horizontal variability of aerosol size 

15 distribution in the vicinity of the SGP site. The Twin Otter aircraft is operated by the 

16 Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) based at 

17 Marina, California [Bluth et al., 1996], and performed 16 research flights out of Ponca 

18 City airport, Oklahoma (32 km east of the SGP central facility). Most of the flight 

19 patterns were centered on the central facility at the SGP site. The details of the Scanning 

20 Mobility Particle Sizer onboard the Twin Otter is given in Wang et al., [2003]. The 

21 major components of the airborne SMPS system are a cylindrical Differential Mobility 

22 Analyzer (TSI Inc., model 3081) and a condensation particle counter (TSI Inc., model 

23 3010). Prior to measurements, the relative humidity (RH) of the aerosol sample was 
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reduced to below 25% inside a Nafion drier. The second SMPS was deployed at the 

2 surface GIF, which is located at the SGP site central facility (36.60°N, 97.48°W, 319 m). 

3 The SMPS at GIF provided continuous measurements of sub-micrometer aerosol size 

4 distributions from May 6 to May 27, except for a few brief periods of instrument 

5 downtime. The only significant difference between the two SMPS systems is that no 

6 Nafion drier was installed in the SMPS at the GIF. The relative humidity (RH) inside the 

7 SMPS at the GIF was lower than the ambient RH as a result of the higher temperature 

8 inside the instrument. The SMPS onboard the Twin Otter measured aerosol size 

9 distribution ranging from 20 nm to 800 nm every 73 seconds; the SMPS at the GIF had a 

10 time resolution of 2 minutes for particles ranging from 25 nm to 1 pm. Both SMPS 

11 systems were carefully calibrated during the IOP using Polystyrene Latex standards. 

12 Data from both SMPS systems were analyzed using the data inversion procedure 

13 described by Collins et al. [2002], A comparison between the size distributions 

14 measured by the two SMPS systems during the IOP is shown in Fig 1. The aerosol size 

15 distributions were averaged from local time 09:57 to 10:12 on May 20, 2003 (local 

16 daylight saving time, GMT-5hrs will be used throughout this paper). During the period, 

17 the Twin Otter flew an overpass above the GIF at an altitude of 200 meters above ground 

18 level (AGL). Two vertical profiles of potential temperature were measured onboard 

19 Balloon-Borne Sounding Systems (BBSS) at 6:30 and 12:30, respectively. The potential 

20 temperature profdes suggest vertical homogeneity from the surface to the altitude of 200 

21 meter (AGL). The relative humidity (RH) inside both instruments was below 25%. The 

22 aerosol size distributions were bimodal, and good agreement was found between the 

23 measurements from the two SMPS. 
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1 2.1. Related aerosol measurements. 

2 Besides the aerosol size distributions, aerosol optical properties measured at the 

3 GIF and onboard the Twin Otter are also used in this study. As part of the long-term 

4 aerosol measurements at the SGP site (since 1996), aerosol scattering, absorption, and 

5 scattering /(RH) have been measured continuously at the Aerosol Observation System 

6 (AOS) trailer, which is located approximately 150 meters from the GIF at the SGP site 

7 [ Sheridan et ah, 2000]. The /(RH) is defined as: 

8 /(RH) 

°/(RH=85%) l®sp( RH=40%) (1) 

9 where <7 (RH=g5%) and cr (RH=40%) are aerosol scattering coefficients at RH of 85% and 

10 40%, respectively. During the IOP, an additional nephelometer and Particle/Soot 

1 1 Absorption Photometer (PSAP) at low RH conditions were set up in the GIF. The details 

12 of the aerosol optical property measurements at the AOS and the GIF during the IOP are 

13 given in Andrews et a!. [2005]. Onboard the Twin Otter aircraft, aerosol scattering, 

14 absorption, and scattering /(RH) were measured by a 3-wavelength nephelometer, a 

15 PSAP, and a humidigraph, respectively. The description of the airborne aerosol optical 

16 property measurements can be found in Schmid et a/. [2005]. The size cut for the 

17 nephelometer and PSAP measurements at the AOS trailer and the GIF were alternated 

18 between 1 pm and 10 pm every 6 minutes. All in-situ instruments onboard the Twin 

19 Otter sampled from a common inlet. The particle transmission efficiency of the inlet is 

20 nearly 100% for particles smaller than 3.5 pm, and 60% for particle of 5.5 pm in diameter 

2 1 for typical Twin Otter velocity of 50 m/s [ Hegg et al., 2005]. 

22 
23 
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2.2. Related trace gas measurements. 

Mixing ratios of trace gases, including CO and S0 2 used in this study, were 
monitored at a NASA mobility facility (http ://smart-commit. gsfc .nasa. gov) located 
approximately 50 meters from the GIF. The CO mixing ratio was measured by a CO 
analyzer (Thermo Environmental Inc, model 48C) that was calibrated with a NIST- 
ratioed CO standard. Data were logged every 5 minutes. The mixing ratio of S0 2 was 
monitored using a S0 2 analyzer (Thermo Environmental Inc, model 43C) that was 
calibrated with a NIST-ratioed S0 2 standard. No chemical zero was installed prior to 
May 18, 2003. After that date, a K 2 C03-impregnated fdter that removes all S0 2 from the 
gas stream was installed and the S0 2 analyzer was zeroed periodically by switching the 
inlet flow to pass through the fdter. As a result, the S0 2 mixing ratios measured after 
May 18 are considered more accurate than those measured before May 18. The 
uncertainties associated with the S0 2 measurements before and after May 18 are 
estimated to be 0.5 ppb and 0.05 ppb, respectively. Because of the large uncertainties in 
the S0 2 measurements relative to the observed mixing ratios, the S0 2 mixing ratios 
measured during the IOP are only used qualitatively in this study. 

3. Results and Discussion 

3. 1. Size distribution measured at surface 

Aerosol number size distributions, aerosol scattering coefficients (Dp<10 pm) at 
three wavelengths, and total particle number concentrations measured at the GIF during 
the IOP are shown in Fig 2. The aerosol size distributions were adjusted to dry 
conditions using the measured hygroscopicity [ Gasparini et al., 2005]. The aerosol 
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1 scattering coefficients were measured at RH mostly below 50% during the IOP. The 

2 precipitation, which is routinely measured at the SGP site, is also presented. During the 

3 IOP, substantial variations were observed in aerosol properties, including aerosol size, 

4 concentration, and scattering coefficients. The aerosol scattering coefficient at 550 nm 

5 ranged from 2.1 to 413 Mm' 1 , with the highest level of 413 Mm" 1 observed on May 8. 

6 This level substantially exceeds the 95 percentile (less than 150 Mm' 1 at 550 nm) of a 

7 four-year record at the SGP site [ Sheridan et al., 2000]. During the same period on May 

8 8, the dry aerosol size distribution exhibited an accumulation-mode diameter of 200 nm, 

9 which was the largest observed at the surface during the IOP. The sources of the aerosol 

10 during this period are discussed later. There were two major precipitation events during 

11 the IOP, May 16 and 24. Following these precipitation events, substantial decreases in 

12 aerosol scattering coefficients were observed; following the precipitation on May 16 and 

13 May 24, the aerosol scattering coefficient at 550 nm decreased from 67 to 9.7 Mm' 1 and 

14 from 62 to 8.1 Mm' 1 , respectively. These decreases are attributed to wet deposition. 

15 High particle number concentrations were frequently observed during the IOP, often 

16 coincident with the observations of nucleation mode particles. 

17 Based on back- trajectory analysis, air masses arriving at the SGP site during the 

18 IOP were grouped into four classes (Figure 2) [ Gasparini et al., 2005]. The back- 

19 trajectories of parcels arriving at the SGP Central Facility 150 m above ground level 

20 (AGL) were calculated using the NOAA HYSPLIT4 model [ Draxler , 1988] at three-hour 

21 intervals for the entire IOP. In order to establish the dominant profdes among these 219 

22 backtrajectories, similar backtrajectories were grouped using the cluster analysis 

23 technique described by Dorling et al. [1992], This analysis resulted in eight distinct back- 
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1 trajectory clusters. Because of the small number of observations associated with several 

2 of the back- trajectory clusters, some of the visually similar clusters were merged into four 

3 classifications that were more statistically-significant. These classes represent 

4 backtrajectories from the Northwest, Midwest, Southeast, and Southwest. Maps showing 

5 the original eight clusters and their grouping into the four back- trajectory classifications 

6 are presented in Fig 3. 

7 The dry sub-micrometer aerosol volume concentration is plotted against the 

8 simultaneously measured CO mixing ratios for each air mass class (Figure 4). The 

9 aerosol volume concentrations are calculated from aerosol number size distributions 

10 measured by the SMPS at the GIF. The aerosol size distributions were adjusted to dry 

1 1 conditions using the hygroscopicity measured by a Tandem Differential Mobility 

12 Analyzer [ Gasparini et al., 2005], and were averaged into 30-minute time intervals. For 

13 air masses from the Southeast, the sub-micrometer aerosol volume concentration was 

14 well correlated with the measured CO mixing ratio, which suggests the aerosols were 

15 mainly produced through combustion processes. During local time 09:12 to 23:15 on 

16 May 8, high sub-micrometer aerosol volume concentration up to 37.2 pm /cm was 

17 observed concurrent with the air mass being from the southeast. During the same period, 

18 the CO mixing ratio reached 560 ppb, and the aerosol scattering coefficient and the 

19 accumulation mode diameter also reached the highest levels observed at the GIF during 

20 the IOP. Back- trajectory analysis indicates the air mass arriving at the SGP during this 

21 period originated in Central America on May 6, at which time extensive wildfires were 

22 reported in Central America. The extensive wildfires are clearly shown in a picture taken 

23 by Moderate Resolution Imaging Spectroradiometer (MODIS) on May 5 th , 2003 
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1 (http://rapidfire.sci.gsfc.nasa.gov/gallery). Using a coupled aerosol, radiation, and 

2 meteorology model, Wang et al. [2006] simulated the transport of the Central American 

3 smoke during April 20 to May 21, 2003. They found the modeled Aerosol Optical 

4 Thickness (AOT) consistently captured the fluctuations of the measured AOT, and the 

5 model-simulated timeline of the smoke events (e.g. May 8-12, 2003) is in good 

6 agreement with those identified from observed AOTs at the SGP site. The simultaneous 

7 observations of high sub-micrometer aerosol volume concentration and CO mixing ratio, 

8 combined with the back trajectory analysis and the modeling results, give strong support 

9 to the attribution of the observation to smoke transported from Central America. Aerosol 

10 compositional analyses by synchrotron x-ray fluorescence of aerosol samples collected at 

11 the surface by a rotating DRUM sampler [ Cahill et al., 1985; Raabe et al., 1988; Cahill 

12 and Wakabayashi, 1993; Bench et al., 2001] showed that this period had high non-soil 

13 potassium indicative of a strong biomass burning signature [ Cahill et al., unpublished 

14 report, 2004], 

15 When the air mass was from the southwest, the sub-micrometer aerosol volume 

16 concentration was also highly correlated with CO mixing ratio. From May 8 23:48 to 

17 May 9 02:12, high CO mixing ratio and sub-micrometer aerosol volume concentration 

18 were again observed. The CO mixing ratio ranged from 345 ppb to 555 ppb, and the sub- 

19 micrometer aerosol volume concentration was 15.8 to 35.1 pm /cm . Except during the 

20 two periods described above, the CO mixing ratio measured during the IOP was mostly 

21 below 350 ppb. 

22 For air masses from the northwest, the CO mixing ratio was mostly between 133 

23 and 200 ppb, and the sub-micrometer aerosol volume concentration varied from 0.55 to 
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1 1.7 jum /cm . The low CO mixing ratio measured is expected as there are no major CO 
sources (such as vehicle emissions from urban areas) towards the northwest, and no lire 
events were reported in the region during the IOP. The large variation in aerosol volume 
concentration and the relatively constant and low CO mixing ratio suggest that the 
aerosol was produced mainly through processes other than combustion. For a small 
fraction of the measurements, the CO mixing ratio reached 300 ppb. Although it is not 
possible to pinpoint the exact reason, the occasional observations of such high CO mixing 
ratio could be due to mixing of the air parcel with vehicle emission from local traffic. 
The Midwest class consisted of more than 50% of the measurements during the IOP. The 
poor correlation between the sub-micrometer aerosol volume concentration and the CO 
mixing ratio is likely due to the following reasons. In Midwest region, CO is mainly 
produced through combustion processes and from cities, which are not significant sources 
of aerosol volume when compared to power plants. Most Midwest power plants emit 
substantial S0 2 but have low CO emission rates. 


3.2. Statistics of aerosol number size distribution during the IOP 

Aerosol size distributions measured at the GIF were fitted with up to three log 
normal functions to derive statistics: 


dN 


1 


N, 


d l°Sio D p ,=i V2^1og 10 cr. 


exp 


(togio^-togio ^) 2 

21 Ogioffi 


( 2 ) 


where N,, D pi and a, are the mode number concentration, geometric mean diameter, and 
geometric standard deviation, respectively. Prior to the fitting, the measured aerosol size 
distributions were adjusted to dry conditions using the measured hygroscopicity 
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[' Gasparini et al, 2005]. The adjusted aerosol size distributions were then averaged into 

2 30-minute intervals. During the IOP, the size distributions observed are described with 

3 three modes: a nucleation mode (D p = 5-30 nm), an Aitken mode (D p =30- 1 00 nm), and an 

4 accumulation mode (D p = 100-500 nm). The statistics of the fitted modes are presented in 

5 Table 1 for each air mass class. As the nucleation modes often exhibited a mode 

6 diameter below 20 nm and were not completely captured by the SMPS, only the 

7 frequency of the appearance of a nucleation mode is considered accurate and is given in 

8 Table 1. The nucleation mode was frequently observed during the IOP for each air mass 

9 class. The mode diameter and the concentration of both Aitken and accumulation modes 

10 showed substantial variations. The average accumulation mode diameters for the 

11 Southeast and Southwest classes are 209 and 199 nm, respectively, which are 

12 substantially larger than those when air masses originated from the Midwest and 

13 Northwest. The larger accumulation mode diameters of the Southeast and Southwest 

14 classes were due to the influence of the extensive fire activities in Central America 

15 described above. Despite the large differences in the accumulation mode diameter and 

16 concentration, the geometric standard deviation of the accumulation mode showed little 

17 variation among the four classes. The average geometric standard deviation of geometric 

18 standard deviation is about 1.50 for all four classes, a value close to the “self-preserving 

19 size distribution” geometric standard deviation of 1.34 [Hinds, 1999]. 

20 

21 3.3. High particle number concentration events observed during IOP 

22 During the IOP, high particle number concentrations, over 15,000 cm' , were 

23 observed 7 times. Figure 5 shows an example of the observed high particle number 
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3 

concentrations, which reached 40,000 cm" at 13:00 on May 7. At the same time, a peak 

2 in SO 2 mixing ratio was recorded, and aerosol number size distribution showed a peak 

3 diameter near or below the 20 mn of the lower size range of the SMPS. The high 

4 concentration of small particles suggests new particle formation through homogeneous 

5 nucleation. For each of the 7 cases, a peak in SO? mixing ratio was observed along with 

6 high particle number concentrations and small peak diameters. Six of the seven cases 

7 were observed during local time from 11:30 am to 15:30, with the remaining one 

8 observed around 17:45. In contrast, there were two periods during which the S0 2 mixing 

9 ratios peaked at over 1 ppb while the total particle number concentrations remained 

10 below 8000 cm" 3 . One of the two periods was at about midnight and the other period was 

11 around 20:00 in local time. As the high particle concentrations were observed mostly 

12 during day time and were coincident with enhanced S0 2 mixing ratios, we think new 

13 particles were generated through homogenous nucleation involving sulfuric acid. The 

14 sulfuric acid vapor is produced through oxidation of S0 2 by OH radical, which is 

15 produced in sunlight. The newly formed particles then grew into the measurement range 

16 of the SMPS through coagulation and condensation. The possible nucleation 

17 mechanisms include binary H 2 S0 4 -H 2 0 and ternary nucleation of H 2 S0 4 , NH 3 , and H?0, 

18 as suggested by previous field studies [ Clarke et al., 1999; Weber et al., 2001; Weber et 

19 al. , 200 1 ; Kulmala et al. , 2002] . 

20 Figure 6 shows the maximum S0 2 mixing ratio and total particle number 

21 concentrations measured during the IOP as a function of the wind direction. High S0 2 

22 mixing ratio and particle number concentrations were observed when the wind was from 

23 the east, south, and southwest. The similar distributions of the S0 2 and particle number 
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concentrations as functions of wind direction support our conclusion that the high particle 

2 concentrations were results of particle formation through homogeneous nucleation 

3 involving sulfuric acid. Several power plants are located to the east, south, and southwest 

4 of the SGP site. To the east of the SGP site are three power plants named Conoco, 

5 Ponca, and Sooner. The Sooner power plant has the highest S02 emission rate of 24 

6 kiloTons/year (kT/year, 1999 EPA emission inventory). All three power plants are 

7 within 45 kilometers of the SGP site. To the Southwest near Enid, Oklahoma is the 

8 Kremlin plant with an emission rate of 6.5kT/yr, and to the south of the SGP site is the 

9 Mustang plant in Oklahoma City, which has a S02 emission rate of 5.0 kT/yr. The 

10 observed excursions of S02 mixing ratio are thus due to the emissions from these or 

1 1 other power plants in the vicinity of the SGP site. When the wind is from east to east- 

12 southeast, the emission from Sooner is likely to dominate the S02 concentration. 

13 

14 3.4 Vertical profiles of size distribution measured onboard the Twin Otter. 

15 In addition to continuous measurements at the surface, a SMPS was deployed on 

16 board the Twin Otter aircraft to characterize the vertical profiles of aerosol size 

17 distribution. During the IOP, 16 research flights were flown out of the Ponca City 

18 airport, Oklahoma. Most of the research flights started between 09:00 to 11:00 in the 

19 morning and ended in the early afternoon. As is typical for a continental site in late 

20 spring, the boundary layer was often stable in the morning, as evidenced by an increasing 

21 potential temperature with increasing altitude. As the surface heating increased, the 

22 boundary layer gradually became well mixed later in the day. Vertical profile 

23 measurements often exhibited higher aerosol number concentration and light scattering 
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1 coefficients within the boundary layer than those inside free troposphere [Andrews et ah, 

2 2004], and rapid transitions were often observed at the top of boundary layer. However 

3 during three flights, unusually high aerosol scattering coefficients were observed at 

4 altitudes over 3000 meters. The aerosol scattering coefficients in these elevated layers 

5 reached near 500 Mm' 'at 530 nm (adjusted to STP condition, 0°C and latm), which is 

6 even higher than those observed inside the Central American plumes at the GIF on May 

7 8. 

8 Figure 7 shows measurements taken during a gradual ascent from 15:40 to 16:12 

9 on May 25. The height of the boundary layer was around 1400 m; a constant potential 

10 temperature within the boundary layer suggests the boundaiy layer was well mixed. In 

1 1 the free troposphere, the potential temperature increased monotonically with increasing 

12 altitude to 5000 m, indicating a stratified free troposphere. The aerosol properties, 

13 including size distribution, scattering and absorption coefficients, are fairly uniform 

14 within the boundary layer. The aerosol number size distribution is unimodal with the 

15 mode diameter around 140 nm and the total particle concentration around 2000 cm' 3 . 

16 The aerosol number concentration and scattering coefficients decreased rapidly from the 

17 top of the boundary layer into the free troposphere, as was observed during most flights 

18 of the IOP. However unlike most of the other flights, two aerosol layers with high 

19 aerosol scattering and absorption coefficients were observed in the free troposphere, at 

20 altitudes of 3150 and 3700 meters. At altitude 3150 m, the aerosol scattering coefficient 

21 at 530 nm exceeded 350 Mm' 1 (adjusted to STP), which is more than 6 times that in the 

22 boundary layer. In contrast, the total particle number concentrations within each of the 

23 two layers were much lower than those in the boundary layer, about 2000 cm' (also 
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1 adjusted to STP). The aerosol size distributions showed a fairly rich structure in the free 

2 troposphere. In the two layers with elevated light scattering coefficients, the 

3 accumulation mode diameters were nearly 400 nm, which accounts for the observed high 

4 aerosol scattering coefficients despite the low number concentrations. 

5 To understand the origin of the elevated aerosol layers, we carried out back- 

6 trajectory calculations using NOAA HYSPLIT4 model. The back-trajectoiy calculation 

7 (Figure 8) indicates the aerosol layers observed originated 1 1 days earlier from Northeast 

8 Asia near the border of China and Russia [Damoah et al., 2004; Jaffe et al., 2004]. 

9 Image 1 is a picture taken by Moderate Resolution Imaging Spectroradiometer (MODIS) 

10 onboard the Terra Satellite on May 14, which indicates extensive forest fires in Russia 

1 1 near the Chinese border. The calculated trajectories over Northeast Asia are at an altitude 

12 of 9 km, which is substantially higher than typical particle injection heights. Despite the 

13 discrepancy in the trajectory height and substantial uncertainties associated with back- 

14 trajectory calculations over such an extended period, the high aerosol scattering 

15 coefficient and large accumulation-mode diameter at high altitude suggest a forest fire as 

16 the only plausible source. The large accumulation-mode diameter, which is up to 400 

17 nm, is likely a result of particle coagulation and condensation of gas phase species during 

18 the transport from northeast Asia to above the SGP site, and is consistent with other 

19 observations of long range transport of aged smoke [ Wandinger et al., 2002; Muryama et 

20 al., 2004; Mattis et al., 2004], 

2 1 Elevated aerosol layers with high scattering and absorption coefficients were also 

22 observed during two other flights on May 27 and 28, 2003. Based on similar 

23 considerations, these elevated layers are also attributed to plumes from the forest fires in 
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the same area in East Asia. The aerosol optical properties averaged over the elevated 

2 layers during the three flights and during the observation of the smoke plume from 

3 Central America at GIF are listed in Table 2. The average aerosol number size 

4 distributions are shown in Figure 9. 

5 The aerosol scattering coefficients observed in the plumes generated by fire 

6 activities ranged from 235 to nearly 500 Mm' 1 . For the elevated plumes from East Asia, 

7 the aerosol size distributions showed large accumulation-mode diameters, near 400 nm. 

8 In contrast, the aerosol observed in plumes originating from wildfire at Central America 

9 exhibited a substantially smaller accumulation-mode diameter, around 200 nm. The 

10 difference in particle sizes observed within the Central American and Asian plumes is 

1 1 possibly due to different fire types and strengths. The difference in the initial scavenging 

12 in pyrocumulus clouds and the scavenging during the long range transport could also lead 

13 to the observed differences in aerosol size distributions. As expected, the aerosols within 

14 the elevated plumes had similar intensive properties: the aerosol single scattering albedo 

15 ranged from 0.96 to 0.97 and the scattering /(RH) ranged from 1.15 to 1.20. The 

16 scattering /(RH) is defined as the ratio of aerosol scattering coefficient at 85% RH to that 

17 at 40% RH. The aerosol single scattering albedo derived from in situ measurements of 

1 8 aerosol scattering and absorption agree with values retrieved from Fidar measurement at 

19 the SGP for the elevated layers [Ferrare et al., 2005b]. The high aerosol single scattering 

20 albedo indicates the aerosol within the layers were relatively non-absorbing. Similar 

2 1 values of single scattering albedo for this smoke were also found using multiwavelength 

22 lidar retrievals for lidar data acquired over Tokyo, Japan and Feipzig, Germany 

23 [Muryama et al., 2004; Mattis et al., 2004]. For plumes from Central America, the 
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1 aerosol single scattering albedo showed a similar value of 0.94, and the scattering /(RH) 

2 measured by AOS was 1.30. 

3 

4 3.5. Study of aerosol evolution through growth law analysis 

5 When the primary aerosol emission and coagulation are negligible, and the 

6 growth of the aerosol particles is controlled mainly by gas-to-particle conversion, the 

7 particle growth rate can be determined from the time series of the aerosol size distribution 

8 [McMurry et al., 1981; Wilson and McMurry, 1981; McMurry and Wilson, 1982], This 

9 technique, which is often referred to as growth law analysis, has been previously applied 

10 to observations of particle growth in power plant plumes, and in ground based and 

1 1 aircraft based studies [McMurry et al., 1981; Wilson and McMurry, 1981; McMurry and 

12 Wilson, 1982; Brock et al., 2002]. The goal of the analysis is to determine the 

13 relationship between the growth rate of particle diameter and the diameter of the particle. 

14 By comparing this relationship to different gas-to-particle conversion mechanisms, it is 

15 possible to determine the mechanism responsible for the observed growth [McMurry et 

16 al., 1981; Wilson and McMurry, 1981; McMurry and Wilson, 1982; Brock et al., 2002]. 

17 The particle diameter growth rate is determined from the rate of the size distribution 

18 change as [McMurry and Wilson, 1982]: 

dD, N'(D t , h )-N'(D pA ) 

dt ! V 'll ", ( n r ) 

20 where D p is the particle diameter, N g the number of particle per cm of air with diameter 

21 greater than D p , and T and tj are the times when measurements were made. n av ( D p ) is 

22 the average size distribution from tj to ty. 
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A discussion of the diameter growth rates under different mechanisms can be found in 
McMurry and Wilson [1982]. For condensation controlled growth, the particle diameter 
growth rate is described by following semi-empirical relation [Bademosi and Liu, 1971]: 


dD„ 


1 


■ oc ■ 


dt D 


1.51 + 2.08 Kn 


p V 


1 .5 1 + Kn(l .08 + 6.04/?) + S.32 j3Kn 2 


(5) 


where Kn = , and J3 = — , a. is the mean free path, D the diffusivity, and c the 

/ U P c / 1 

mean thermal velocity of the condensing species. For condensation controlled growth, 
the diameter growth rate is independent of particle size in the free molecular regime 
( D p 2 ), and decreases proportional to 1 ID in the continuum regime (D » /t). As 

the SGP site was in the vicinity of several power plants, and high SO 2 mixing ratio was 
frequently observed during the IOP, one of the likely condensing species is the H 2 SO 4 
formed through the gas phase oxidation of SO 2 . The diffusivity of hydrated sulfuric acid 
vapor in air has been experimentally determined to be 0.077 cm's' [Hanson and Eisele, 
2000 ], 

If the particle growth is due to volume controlled reaction/growth (i.e. aqueous 
reaction within hydrated particles, gas/particle partitioning), the corresponding growth 
rate is expressed as [McMurry and Wilson, 1982]: 


dl\ 

dt 


ccD p 


( 6 ) 


Using the continuous aerosol size distribution measured at the GIF, growth law analyses 
were performed in periods during which the assumptions for growth analysis were 
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1 


satisfied. The assumptions are (1) during the chosen periods, the primary particle 

2 emissions and coagulation were not significant, as evidenced by size distributions and a 

3 constant total particle number concentration; (2) the growth of the particles over the 

4 measurement period was substantial, such that the growth rates derived are statistically 

5 significant; (3) there is no significant dilution or change of air mass, as evidenced by a 

6 constant total particle number concentration. Based on these criteria, four cases are 

7 chosen for growth law analysis; the results are presented below. 

8 

9 3.5.1. Case I of growth law analysis, May 16, 2003. 

10 The first case analyzed is from local time 14:15 to 21:15 on May 16. During the 

1 1 period, the sky was overcast with occasional light drizzle, and the aerosol size 

12 distribution showed continuous particle growth (Figure 10). The cloud was dominated by 

13 mid-level convective precipitating systems with some low-level stratus. The cloud 

14 coverage was 100%, and the base height of the stratus cloud was around 300 meters 

15 (AGL). Based on the measured precipitation rate, the calculated scavenging coefficient is 

16 less than 9e-3/hr for particles ranging from 20nm to 1pm in diameter, and the effect of 

17 scavenging from the drizzle is negligible on particle size distribution during this period. 

18 The aerosol total number concentration, which was measured by a condensation particle 

19 counter (CPC), remained constant at around 2000 cm' , suggesting primary particle 

20 emission, coagulation, and wet scavenging are negligible. Figure 1 1 presents the particle 

21 size distributions at the beginning and the end of the period; size distributions were 

22 bimodal throughout the period. Substantial growth was observed for both modes. The 

23 mode diameter of the Aitken mode grew from 54 to 76 nm and the mode diameter of the 
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accumulation mode increased from 127 to 208 nm. Figure 11 also shows the particle 
growth rates derived using Eqn. 3. The uncertainties in derived growth rate are estimated 
based on counting statistics and sizing accuracy of the SMPS system, and the variation of 
particle concentration derived from the total particle concentrations. A least squares fit of 
the diameter growth rate shows the size dependence of the growth rate is well represented 
by a combination of condensation growth and volume controlled reaction: 


dD 

dD 


dD 


p 

_ p 

+ p 


dt 

dt 

C 

dt 

V.R. 

dD p 

l n m 

1 I 

f 

1.51 + 2.08 Kn 

dt 

c D p \ 

yl.51 + Kn( 1 .08 + 6.04/?) + 8.32 /3Kn J 

dD,, 






i =b-D 



dt 

V.R. 

P 



( 7 ) 


where p 


dt 


and dDp 


dt 


are the fitted particle diameter growth rates due to 


condensation and volume controlled reactions, respectively, and a and b are fitted 
coefficients. The fractional contribution of the volume controlled reaction to the increase 
of the aerosol volume concentration is given by: 


fv.R. = ' 


r'% D 2 dD P 

■’Ai p dt 

n cn 

V.R. 

,[D p ,t)dD p 

\ Dpl D 2 

J D pl P 

1 

& 

+ dDp 


n ar( D p ,t)dD p 

dt 

V.R. dt 

cj 


( 8 ) 


where n av ( D p , t j is the average particle size distribution during the period. For this case, 


the fv.R. is 94%, suggesting the increase of aerosol volume concentration was dominated 
by volume controlled reactions. 
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1 


Size resolved aerosol chemical composition was inferred from aerosol 

2 hygroscopicity measured by a Tandem Differential Mobility Analyzer system at the GIF 

3 during the IOP [ Gasparini et al., 2005; Gasparini et al., 2004], In their analysis, 

4 Gasparini et al. assumed that the particles consisted of ammonium sulfate and a group of 

5 low-hygroscopic species, which showed hygroscopic growth factors (ratio of particle 

6 diameter at 85% RH to that at 40%) less than 1.1 [ Gasparini et al., 2004], As the 

7 hygroscopic growth of ammonium sulfate is well known, size resolved sulfate volume 

8 fraction can be derived from the measured aerosol hygroscopicity [ Gasparini et al., 

9 2004]. The assumption of sulfate being the only inorganic hygroscopic species is 

10 supported by: (1) the contribution of sea salt to the aerosol mass is negligible given the 

1 1 inland location of the SGP site; (2) the chemical composition of sub-micrometer aerosol 

12 measured by a Particle Into Liquid Sampler (PILS) at the GIF showed the contributions 

13 of both sea salt and nitrate were small, and (3) the sulfate was often completely 

14 neutralized by ammonium. Based on the sulfate fraction derived from aerosol 

15 hygroscopicity, at the beginning of the first case, 14:15 on May 16, the volume 

16 concentration of sulfate and low-hygroscopic species are calculated as 1.01 and 0.64 

17 pm /cm , respectively. The volume concentrations of sulfate and low-hygroscopic 

18 species increased to 6.76 and 2.92 pm /cm at 21:15. Sulfate represents 72% of the total 

19 increase in aerosol volume concentration. 

20 The dominance of sulfate in the increase of aerosol volume concentration, 

21 together with the volume controlled reaction mechanism inferred from growth law 

22 analysis, suggests aqueous oxidations of S0 2 by H 2 0 2 and/or 0 3 within hydrated aerosol 

23 particles a possible source of the observed sulfate. The conclusion is also supported by 
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1 the high RH observed during this period, which is averaged at 94% at surface. As a 

2 result of adiabatic cooling, the ambient RH increased with increasing altitude, and 

3 reached 100% at the cloud base, indicating enhanced SO 2 oxidation rate due to higher 

4 liquid water content in hydrated particles. The nearly 100% cloud coverage and the 

5 bimodal size distributions observed during the period also suggest possible in cloud 

6 production of sulfate through aqueous oxidation of SO 2 . As suggested by Hoppel et al., 

7 [1986], the larger cloud condensation nuclei (CCN) grew through aqueous phase 

8 reactions in clouds to form the larger accumulation mode, while the small non-CCN grew 

9 at slower pace through condensation process, and became the Aitken mode. 

10 Measurements onboard the BBSS at 12:30 and 18:30 showed uniform potential 

1 1 temperature and water vapor mixing ratio within the boundary layer, indicating a well 

12 mixed boundaiy layer. After particles grew through aqueous reaction in high RH region 

13 and in clouds, they were transported down to the surface through detrainment and mixing 

14 process, which explains the observed particle growth at the GIF. For in cloud production 

15 of sulfate dominated by the aqueous reaction of S0 2 and H 2 O 2 , the particle volume 

16 growth rate is proportional to the product of the cloud droplet volume and the fraction of 

17 time during which particles are activated into cloud droplets: 


dV _ n d2 dD p 
dt 2 p dt 


6 


(9) 


19 where A/ is the droplet diameter after particles are activated and f is the fraction of time 

20 when particles are activated. The particle diameter growth rate is derived as: 



( 10 ) 
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1 Based on condensation of supersaturated water vapor, adiabatic cloud models often 

2 predict narrow droplet size distributions in clouds (i.e. CCN of different diameters are 

3 activated and form cloud droplets with similar diameters). However, many field studies 

4 have shown cloud droplet size distributions observed had higher dispersions than those 

5 predicted by the models [. Martin et al., 1994; Hudson and Svensson 1995; Politovich 

6 1993; Hudson and Yum 1997]. Larger CCN are likely activated and form larger cloud 

7 droplets. Moreover, as aerosol particles experience different water vapor 

8 supersaturations due to inhomogeneity of updraft velocity and turbulence, larger particles 

9 were more likely to be activated into cloud droplets and the fraction of time (f ) during 

10 which they remain as cloud droplets is higher than that of smaller particles. Both larger 

1 1 droplet size and higher f lead to higher aqueous sulfate production rate, and possibly 

12 result in increased particle diameter growth rate as increasing diameter, as shown by 

13 growth law analysis. Similar particle growth law showing combination of condensational 

14 growth and volume controlled reactions have been observed in previous field studies 

15 [McMurry et al., 1981; McMurry and Wilson, 1982] including a study of aerosol growth 

16 during a winter fog event [Ulevicius et al., 1994], 

17 Besides sulfate, the volume concentration of low-hygroscopic species also 

18 increased substantially from 0.64 to 2.92 pm /cm . As the contributions from mineral 

19 dust and elemental carbon to sub-micrometer aerosol volume concentration are small at 

20 the SGP, we think most of the low-hygroscopic species were organics, which were likely 

21 formed through partitioning of secondary organics into the particle phase. As the volume 

22 growth resulting from the partitioning of secondary organics increases with the increasing 

23 volume of pre-existing organic species within the particle [Pankow 1994], the rate of 
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1 particle growth due to gas/particle partitioning has the same size dependence (growth 

2 law) as those under volume controlled reactions. Size resolved chemical speciation, 

3 which were not available during the IOP, would be able to provide more information 

4 regarding the mechanisms responsible for the increase of low-hygroscopic volume 

5 concentration. 

6 

7 3.5.2. Case II of growth law analysis, May 17 th . 

8 The second case analyzed is from 4:45 to 12:45 on May 17. The sky was overcast 

9 with stratus cloud base height at around 500 meters (AGL). The average RH reached 

10 86% at the surface during this period. The aerosol size distribution, sub-micrometer 

11 volume concentration, and total number concentration are presented in Figure 10. 

12 Similar to the first case, aerosol size distributions were bimodal. Figure 12 shows size 

13 distributions at the beginning and the end of the period, along with the derived diameter 

14 growth rate. The particle growth rate is again well represented by a combination of 

15 condensation growth and volume controlled reactions. While the volume controlled 

16 reactions and/or gas/particle partitioning still dominated the overall particle growth, their 

17 contributions to the total increases in aerosol volume concentration decreased from 94% 

18 in the first case to 64%. Based on aerosol hygroscopicity measurements, the sulfate 

19 volume concentration increased 4.2 pm /cm from 5.6 pm /cm to 9.8 pm /cm , which 

20 represents 71% of the total increase in volume concentration. Given the high RH and 

21 nearly 100% cloud coverage, and the dominance of sulfate in the increased aerosol 

22 volume concentration, the sulfate, at least partially, was produced through aqueous 

23 oxidation of SO 2 . Without further information, we can not determine the species that 
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1 were involved in the condensational growth. One of the possible condensing species is 

2 H2SO4, produced through gas phase reaction of SO2 and OH. 

3 No significant particle growth was observed and the sub-micrometer aerosol 

4 volume concentration remained nearly constant before 07:15. Vertical profiles of 

5 potential temperature and water vapor mixing ratio showed a stable boundary layer at 

6 06:30. The boundary layer became well mixed at 12:55, which is evidenced by the 

7 uniform potential temperature and water vapor mixing ratio within the boundary layer. 

8 The constant sub-micrometer aerosol volume concentration observed could be explained 

9 by the following. First, aqueous production of sulfate is expected to be most significant 

10 within the clouds and in the regions below clouds where RH is nearly 100%, and a stable 

1 1 boundary layer prevented aerosol processed by the clouds and within high RH regions, if 

12 any, to reach the surface level; second, aqueous production of sulfate was most likely 

13 dominated by the reaction of S0 2 and H 2 O 2 . As the H 2 O 2 is produced in sunlight, we do 

14 not expect substantial aqueous production of sulfate in the early morning. After 7:15, the 

15 boundary layer gradually became well mixed due to the solar heating of the surface. The 

16 aerosol particles processed in clouds and regions with high RH were then transported 

17 down to the surface, which explains the observed particle growth at the GIF. 

18 

19 3.5.3. Case III and Case IV 

20 Two additional cases were analyzed for particle diameter growth rates. Both 

21 cases were in the morning with cloud coverage over 90%. Case III is from 7:45 to 10:45 

22 on May 20, and Case IV is from 8:15 to 11:15 on May 25. The size distributions and 

23 calculated diameter growth rates are given in Figures 13 and 14. Due to the shorter 
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1 growth time, the particle growth observed during cases III and IV is substantially less 

2 than that during cases I and II. As a result, the uncertainties in derived diameter growth 

3 rates are significantly higher. Least square fits of the diameter growth rates indicate the 

4 particle growth in case III was due to a combination of condensation and volume 

5 controlled reactions, whereas the particle growth in case IV is a result of volume 

6 controlled reactions only. For case IV, the uncertainty in derived diameter growth rate 

7 increases significantly with decreasing particle diameter, especially for diameters less 

8 than 100 nm. As the diameter growth rate is the highest at the small particle size range 

9 under condensation processes, we can not exclude the possibility of condensational 

10 growth because of the large uncertainties in diameter growth rates at small particle sizes. 

1 1 However, for both cases, the uncertainties in diameter growth rates are small for particles 

12 larger than 100 nm, and the volume controlled reactions dominate the particle growth, 

13 especially for particles having diameters larger than 100 nm. 

14 The results of growth law analyses are summarized in Table 3. For all four cases 

15 analyzed, the volume controlled reactions dominate the overall increases in aerosol 

16 volume concentration. Based on chemical composition inferred from hygroscopicity 

17 measurements, ammonium sulfate contributes to 50%-72% of the total aerosol volume 

18 concentration increase. Extensive cloud coverage and high RH below the clouds suggest 

19 the sulfate was likely produced through aqueous oxidation of S0 2 in cloud droplets and 

20 hydrated aerosol particles. The produced sulfate was then transported down to the 

21 surface as a result of mixing process within the boundaiy layer. Besides sulfate, low- 

22 hygroscopic species also contributed significantly to the increases of aerosol volume 

23 concentration (28%-50%). Based on the particle growth law, the increases in volume 
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1 


concentration of low-hygroscopic species were likely due to the partitioning of the 

2 secondary organic species into the particle phase. For all four cases analyzed, the 

3 condensational growth had small contributions to the overall particle growth. Given the 

4 high cloud coverage during the periods analyzed, it is expected that the growth due to 

5 condensation was less important as the actinic flux and photochemical reaction rates were 

6 low. For comparison, it would be interesting to study the particle growth rates when 

7 cloud coverage was low. Unfortunately, no such case satisfying the assumptions for 

8 growth law analysis was available during the IOP. Future work will include studies of 

9 particle growth rates and mechanisms under different meteorological conditions and at 

10 various locations. 

11 

12 4. Conclusions 

13 Aerosol size distributions were measured continuously at the GIF and onboard the 

14 Twin Otter aircraft at the SGP site during May 2003, as part of the Aerosol Intensive 

15 operation period of the DOE ARM program. Large variations in aerosol properties, 

16 including size distribution and optical properties were observed. Using simultaneous 

17 trace gas measurements and back trajectory analyses, the aerosol properties were found to 

18 be influenced by both local and long-range sources. The aerosol from the southeast was 

19 found to be influenced by wild fires in Central America. Aerosol with high scattering 

20 coefficients (up to 500 Mm' 1 ), which can be traced back to forest fires in East Asia, was 

21 also observed at altitudes above 3000 meters. These aerosols had large mode diameters 

22 near 400 nm in diameter, and were most likely aged in plumes during their long transport 

23 over the Pacific Ocean from East Asia. High particle concentrations observed during the 
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1 IOP were well correlated with high SO 2 mixing ratios, suggesting nucleation involving 

2 sulfuric acid is likely the main source of newly formed particles. The high SO 2 

3 concentrations were observed when the wind was from the east, south, and southeast, 

4 where several power plants are located. 

5 Detailed growth law analyses were carried out for four periods, during which 

6 nearly 100% cloud coverage was observed. The results indicate the particle growth in the 

7 four cases was due to a combination of condensation and volume controlled reactions, 

8 and the volume controlled reactions dominated the overall growth. Besides sulfate, low- 

9 hygroscopic species, likely secondary organics, also contributed substantially to the 

10 particle growth. For two cases, the low-hygroscopic species accounted for nearly 50% 

1 1 percent of the total increases in aerosol volume concentration. The growth laws and 

12 meteorological conditions indicate the sulfate was mainly produced through aqueous 

13 oxidation of SO 2 in cloud droplets and hydrated aerosol particles, and the increases in 

14 low-hygroscopic volume concentration were likely due to partitioning of secondary 

15 organic species into the particles. The results from growth law analyses suggest when the 

16 cloud coverage is high, the fractional contribution of condensational growth to overall 

17 particle growth are small at the SGP. 


29 



5. Acknowledgments 

This paper has been authored with support from the Office of Biological and 
Environmental Research under contract number DE-AC02-98CH 10866 and DE-AC06- 
76RL01830 (ITF#5788) with the U.S Department of Energy. We thank Steve Schwartz, 
Larry Kleinman, Yin-Nan Lee, and Peter Damn for their comments on this manuscript, 
and Mike Jensen for provide meteorological information during the IOP. The authors 
also acknowledge the US Dept, of Energy ARM program for providing support at the 
SGP site and for aircraft measurements. 


30 



Reference: 


Andrews, E., P. J. Sheridan, J. A. Ogren, and R. Ferrare (2004), In situ aerosol profiles 
over the Southern Great Plains cloud and radiation test bed site: 1. Aerosol optical 
properties, J. Geophys. Res., 109( D6). Art. No. D06208. 

Andrews, E., et al. (2005), Comparison of methods for deriving aerosol asymmetry 
parameter,/. Geophys. Res., Accepted. 

Bademosi, F., and B. Y. H. Liu (1971), Diffusion charging of Rnudsen aerosols - I. 
(Experiment) and II. (Theory). Publications 155, 156 pp, Particle Technology Laboratory, 
Minnesota. 

Bench, G., P. G. Grant, D. Ueda, S.S. Cliff, K.D. Perry, and T. A. Cahill (2001), The use 
of STIM and PESA to respectively measure profiles of aerosol mass and hydrogen 
content across Mylar rotating drum impactor samples, Aerosol Sci. Technol., 36:6 42-651. 

Bergin, M. EL, S. E. Schwartz, R. N. Halthore, J. A. Ogren, and D. L. Hlavka (2000), 
Comparison of aerosol optical depth inferred from surface measurements with that 
determined by Sun photometry for cloud-free conditions at a continental US site, J. 
Geophys. Res., 105( D5). 6807-6816. 

Bluth, R. T., P. A. Durkee, J. H. Seinfeld, R. C. Flagan, L. M. Russell, P. A. Crowley, 
and P. Finn (1996), Center for Interdisciplinary Remotely-Piloted Aircraft Studies 
(CIRPAS), B. Am. Meteorol. Soc., 77(1 1). 2691-2699. 

Brock, C. A., et al. (2002), Particle growth in the plumes of coal-fired power plants, J. 
Geophys. Res., 707(D12). Art. No. 4155. 

Cahill, T. A., C. Goodart, J. W. Nelson, R. A. Eldred, J. S. Nasstrom, and P. J. Feeney 
(1985), Design and evaluation of the drum impactor. Proceedings of International 
Symposium on Particulate and Multi-phase Processes. Teoman Ariman and T. Nejat 
Veziroglu, Editors. Hemisphere Publishing Corporation, Washington, D.C. 2:319-325. 

Cahill, T. A. and P. Wakabayashi (1993), Compositional analysis of size-segregated 
aerosol samples. Chapter in the ACS book Measurement Challenges in Atmospheric 
Chemistry. Leonard Newman, Editor. Chapter 7, Pp. 211-228. 

Clarke, A. D., V. N. Kapustin, F. L. Eisele, R. J. Weber, and P. H. McMurry (1999), 
Particle production near marine clouds: Sulfuric acid and predictions from classical 
binary nucleation, Geophysical Research Letters, 26(16). 2425-2428. 

Collins, D. R., R. C. Flagan, and J. H. Seinfeld (2002), Improved inversion of scanning 
DMA data, Aerosol Sci. Technol., 36(1). 1-9. 


31 



Damoah, R., et al. (2004), Around the world in 17 days - hemispheric-scale transport of 
forest fire smoke from Russia in May 2003, Atmospheric Chemistry And Physics, 4 1311- 
1321. 

Dorling, S. R., T. D. Davies, and C. E. Pierce (1992), Cluster- Analysis - A Technique For 
Estimating The Synoptic Meteorological Controls On Air And Precipitation Chemistry - 
Method And Applications, Atmospheric Environment Part A-General Topics, 26(14). 
2575-2581. 

Draxler, R. R. (1988), Hybrid Single-Parcel Lagrangian Integrated Trajectories 
(HYSPLIT): Model Description, NOAA Tech. Memo., Natl. Oceanic and Atmos. 
Admin., Washington, D. C. 

Feingold, G., W. L. Eberhard, D. E. Veron, and M. Previdi (2003), First measurements of 
the Twomey indirect effect using ground-based remote sensors, Geophys. Res. Lett., 
30(6). 

Ferrare, R. A., G. Feingold, S. Ghan, J. A. Ogren, S. E. Schwartz, B. Schmid, and P. 
Sheridan (2005a), Preface to Special Section on the DOE ARM May 2003 Aerosol 
Intensive Operation Period (IOP), J. Geophys. Res., submitted. 

Ferrare, R.A., D.D. Turner, M. Clayton, B. Schmid, J.Redemann, D. Covert, R. Elleman, 
J. Ogren, E. Andrews, J.E.M. Goldsmith, H.H. Jonsson (2005b), Evaluation of Daytime 
Measurements of Aerosols and Water Vapor made by an Operational Raman Lidar over 
the Southern Great Plains, J. Geophys. Res. 2005JD005837, in press, this issue 

Gasparini, R., R. J. Li, and D. R. Collins (2004), Integration of size distributions and size- 
resolved hygroscopicity measured during the Houston Supersite for compositional 
categorization of the aerosol, Atmospheric Environment, 55(20). 3285-3303. 

Gasparini, R., R. J. Li, D. R. Collins, and R. A. Ferrare (2005), Aerosol hygroscopicity 
and volatility measured at the ARM SGP site and their application to investigate 
composition and evolution, J. Geophys. Res., Accepted. 

Hanson, D. R., and F. Eisele (2000), Diffusion of H2S04 in humidified nitrogen: 
Hydrated H2S04, J. Phys. Chem. A, 104(8). 1715-1719. 

Hegg, D. A., D. S. Covert, H. Jonsson, and P. A. Covert (2005), Determination of the 
transmission efficiency of an aircraft aerosol inlet, Aerosol Sci. Technol., Submitted. 

Hinds, W. C. (1999), Aerosol Technology, pp. 260-265 and 285-288 pp., John Wiley & 
Sons, New York. 

Hoppel, W. A., G. M. Frick, and R. E. Larson (1986), Effect Of Nonprecipitating Clouds 
On The Aerosol Size Distribution In The Marine Boundary-Layer, Geophysical Research 
Letters, 13(2). 125-128. 


32 



Hudson, J. G., and G. Svensson (1995), Cloud Microphysical Relationships In California 
Marine Stratus, Journal Of Applied Meteorology, 54(12). 2655-2666. 

Hudson, J. G., and S. S. Yum (1997), Droplet spectral broadening in marine stratus, 
Journal Of The Atmospheric Sciences, 54(22). 2642-2654. 

Intergovernmental Panel on Climate Change (2001), Climate Change 2001: The 
Scientific Basis, Cambridge University Press, New York. 

Jaffe, D., I. Bertschi, L. Jaegle, P. Novelli, J. S. Reid, H. Tanimoto, R. Vingarzan, and D. 
L. Westphal (2004), Long-range transport of Siberian biomass burning emissions and 
impact on surface ozone in western North America. Geophys. Res. Lett., 31, 
doi: 10.1 029/2004GL020093 . 


Kim, B. G., S. E. Schwartz, M. A. Miller, and Q. L. Min (2003), Effective radius of cloud 
droplets by ground-based remote sensing: Relationship to aerosol, J. Geophys. Res., 
108(023). 

Kulmala, M., P. Korhonen, I. Napari, A. Karlsson, H. Berresheim, and C. D. O'Dowd 
(2002), Aerosol formation during PARFORCE: Ternary nucleation of H2S04, NH3, and 
H20, Journal Of Geophysical Research-Atmospheres, 107( D19). Art. No. 8111. 

Martin, G. M., D. W. Johnson, and A. Spice (1994), The Measurement And 
Parameterization Of Effective Radius Of Droplets In Warm Stratocumulus Clouds, 
Journal Of The Atmospheric Sciences, 57(13). 1823-1842. 

Mattis I., A. Ansmann, U. Wandinger, D. Muller (2003), Unexpectedly high aerosol load 
in the free troposphere over central Europe in spring/summer 2003, Geophys. Res. Lett., 
30 (22), 2178, doi:10.1029/2003GL018442. 

Murayama T., D. Muller, K. Wada, A. Shimizu, M. Sekiguchi, T. Tsukamoto (2004), 
Characterization of Asian dust and Siberian smoke with multi-wavelength Raman lidar 
over Tokyo, Japan in spring 2003, Geophys. Res. Lett., 31, L23103, 

doi: 10.1 029/2004GL02 1105. 

Pankow, J. F. (1994), An Absorption-Model Of The Gas Aerosol Partitioning Involved In 
The Formation Of Secondary Organic Aerosol, Atmospheric Environment, 28(2). 189- 
193. 

Politovich, M. K. (1993), A Study Of The Broadening Of Droplet Size Distributions In 
Cumuli , Journal Of The Atmospheric Sciences, 50(14). 2230-2244. 

Raabe, O. G., D. A. Braaten, R. L. Axelbaum, S. V. Teague, and T. A. Cahill (1988), 
Calibration Studies ofthe DRUM Impactor. J. Aerosol Sci., 79(2): 183-1 95. 

Schmid, B., et al. (2005), How well do state-of-the-art techniques measuring the vertical 
profde of tropospheric aerosol extinction compare? J. Geophys. Res., Accepted. 


33 



Sheridan, P. J., D. J. Delene, and J. A. Ogren (2001), Four years of continuous surface 
aerosol measurements from the Department of Energy's Atmospheric Radiation 
Measurement Program Southern Great Plains Cloud and Radiation Testbed site, J. 
Geophys. Res., 706(D18). 20735-20747. 

Turner, D. D., R. A. Ferrare, and L. A. Brasseur (2001), Average aerosol extinction and 
water vapor profiles over the southern great plains, Geophys. Res. Lett., 28(23). 4441- 
4444. 

Turner, D. D., R. A. Ferrare, L. A. H. Brasseur, and W. F. Feltz (2002), Automated 
retrievals of water vapor and aerosol profiles from an operational Raman lidar, J. Atmos. 
Oceanic Technol., 19( 1). 37-50. 

Ulevicius, V., S. Trakumas, and A. Girgzdys (1994), Aerosol-Size Distribution 
Transformation In Fog, Atmospheric Environment, 28(5). 795-800. 

Wandinger et al. (2002), Optical and microphysical characterization of biomass burning 
and industrial pollution aerosols from multiwavelength lidar and aircraft measurements, 
J. Geophys. Res., 107, doi:10.1029/2000JD000202. 


Wang, J., R. C. Flagan, and J. H. Seinfeld (2003), A differential mobility analyzer 
(DMA) system for submicron aerosol measurements at ambient relative humidity, 
Aerosol Sci. Technol., 37(1). 46-52. 

Wang. J., S. A. Christopher, U. S. Nair, J. S. Reid, E. M Prins, J. Szykman, and J. L. 
Hand (2006), Mesoscale modeling of Central American smoke transport to the United 
States: 1. “Top-down” assessment of emission strength and diurnal variation impacts, J. 
Geophys. Res., 777(D05S17), doi: 10.1 029/2005 JD0064 16. 

Weber, R. J., P. H. McMurry, R. L. Mauldin, D. J. Tanner, F. L. Eisele, A. D. Clarke, and 
V. N. Kapustin (1999), New particle formation in the remote troposphere: A comparison 
of observations at various sites, Geophysical Research Letters, 26(2). 307-310. 

Weber, R. J., et al. (2001), Measurements of enhanced H2S04 and 3-4 nm particles near 
a frontal cloud during the First Aerosol Characterization Experiment (ACE 1), Journal Of 
Geophysical Research-Atmospheres, 106( D20). 24107-24117. 


34 



Figure Captions 


Figure 1. Comparison of the aerosol size distributions measured by the SMPS onboard 
the Twin Otter (solid line) to that measured at the GIF (dashed line). The aerosol size 
distributions were averaged from local time 09:57 to 10:12 on May 20, 2003. During this 
period, the Twin Otter flew over the GIF at an altitude of 200 meters above ground level 
(AGL). 

Figure 2. Measurements from May 6 to May 27, 2003 at the SGP surface sites, (a) 
Aerosol number size distributions measured by SMPS, adjusted to diy conditions, (b) 
Aerosol scattering coefficients (D p < 10 pm) measured by a 3 -wavelength nephelometer at 
RH mostly below 50%. (c) Precipitation and total particle number concentration 
measured by a condensation particle counter, (d) Classification of air mass based on 
back- trajectory analyses. Note in (a), the maximum concentration in the color scale is set 
at 20,000 cm" ' in order to show details of the size distribution evolution. 

Figure 3. Eight clusters of 219 back- trajectories and their grouping into four back- 
trajectory classes. The back- trajectories were calculated for air parcels arriving at the 
SGP site 150 m AGL every three hours during the IOP. 

Figure 4. Sub-micrometer aerosol volume concentration plotted again simultaneous 
measurements of CO mixing ratio for each air mass class. 

Figure 5. Measurements during a nucleation event from 07:00 to 19:00 on May 7. (a) 
Dry aerosol number size distributions with the color indicating the number concentration, 
(b) SO? mixing ratio and number concentration of particles having diameters greater than 
10 nm. 

Figure 6. Maximum SO 2 mixing ratio and particle number concentrations measured 
during the IOP as functions of wind direction. 

Figure 7. Vertical profiles of measurements taken during an ascent from 15:40 to 16:12 
(local time) on May 25, 2003. (a) Relative humidity and potential temperature, (b) 
Aerosol scattering and absorption coefficients at 530 nm, adjusted to STP (0°C and 1 
atm), (c) Number concentration of particles having diameters greater than 10 nm, 
measured by a Condensation Particle Counter, concentration adjusted to STP, and (d) 
Aerosol number size distribution adjusted to STP. 

Figure 8. Result of back- trajectory analysis for elevated aerosol layer observed above the 
SGP on May 25, 2003. 

Figure 9. Aerosol number size distributions averaged over elevated layers during three 
flights and during the observation of plumes from Central America at the GIF. 
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Figure 10. Measurements from May 16, 12:00 to May 17, 15:00. (a) Ambient 
temperature and RH. (b) Sub-micrometer aerosol volume concentrations calculated from 
aerosol number size distributions, (c) Aerosol number size distributions, with the color 
indicating the number concentration, (d) Total particle number concentration (D p > 10 nm) 
and precipitation. 

Figure 11. (a) Aerosol size distributions at the beginning and the end of the first case for 
growth law analysis, and (b) derived particle growth rate and the least square fit to 
growth rates of condensational processes (Eqn. 5) and volume based reactions (Eqn. 6). 

Figure 12. Same as figure 11 for the measurements from 04:45 to 12:45, May 17 (the 
second case). 

Figure 13. Same as figure 11 for the measurements from 07:45 to 10:45, May 20 (the 
third case). 

Figure 14. Same as figure 11 for the measurements from 08:15 to 11:15, May 25 (the 
fourth case). 
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Table 1. Statistics of fitted aerosol number size distributions during the IOP. 


Mode 

Parameters 

Midwest 

Southeast 

Southwest 

Northwest 

Nucleation 

Occurrence 

47% 

69% 

35% 

52% 

Aitken 

Occurrence 

83% 

76% 

71% 

92% 

D pi (nm) 

69±29 

90±53 

51±1 1 

62±30 

Nj (cm 3 ) 

2304±1620 

1368±693 

6952±6022 

3369±3358 

<?i 

1.73±0.36 

1.84±0.47 

1.65±0.36 

1.46±0.27 

Accumulation 

Occurrence 

93% 

100% 

71% 

91% 

D p:i (nm) 

165±40 

209±33 

199±13 

151±38 

Ni (cm 3 ) 

1202±913 

2301±1368 

1841±1263 

749±794 

<Ji 

1.52±0.17 

1.47±0.13 

1.51±0.04 

1.53±0.12 


# of size dist. 

481 

51 

17 

205 


Table 2. Aerosol optical properties averaged over elevated layers during three flights and during 
the observation of plumes from Central America at the GIF. 


Date 

Altitude 

(meter) 

_ (i 

Osp 

(Mm 1 ) 

_ d 

^ap 

(Mm 1 ) 

Single 

Scattering Albedo 

/(RH) 

(85%/40%) 

5/25/2003 

3568 

497. 2 a 

21.1 a 

0.96 a 

1.20 c 

5/27/2003 

3041 

235. 0 a 

9.2 a 

0.96 a 

1.1 5 C 

5/28/2003 

4132 

292.T 

9.4 a 

0.97 a 

1.20 c 

5/9/2003 

Ground (SGP) 

359. 8 b 

21.9 b 

0.94 b 

1.30 b 


at 530nm. b at 550nm. c at 540nm. d at STP (0°C and 1 atm= 10 1,325 Pa) 
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Table 3. Summary of particle growth law analysis using size distributions measured at 
the GIF during the IOP. 


Date 

Period 
(local time) 

Air 

mass 

Growth 

mechanisms 

/v.R. 

AV 

(pm 3 /cm 3 ) 

A F(sulf 
ate)d V 

Cloud 

coverage: 

May 16, 2003 
(Case I) 

14:15-21:15 

Midwest 

Condensation & 
Vol. reaction 

94% 

8.03 

72% 

100% 

May 17, 2003 
(Case II) 

4:45-12:45 

Midwest 

Condensation & 
Vol. reaction 

64% 

5.91 

71% 

100% 

May 20, 2003 
(Case III) 

7:45-10:45 

Northwe 

st 

Condensation & 
Vol. reaction 

81% 

1.29 

50% 

90%- 

100% 

May 25, 2003 
(Case IV) 

8:15-11:15 

Midwest 

Vol. reaction 

100% 

1.56 

57% 

90%- 

100% 
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Figure 1. Comparison of the aerosol size distributions measured by the SMPS onboard 
the Twin Otter (solid line) to that measured at the GIF (dashed line). The aerosol size 
distributions were averaged from local time 09:57 to 10:12 on May 20, 2003. During this 
period, the Twin Otter flew over the GIF at an altitude of 200 meters above ground level 
(AGL). 
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Figure 2. Measurements from May 6 to May 27, 2003 at the SGP surface sites, (a) 
Aerosol number size distributions measured by SMPS, adjusted to dry conditions, (b) 
Aerosol scattering coefficients (D p <\0 pm) measured by a 3 -wavelength nephelometer at 
RH mostly below 50%. (c) Precipitation and total particle number concentration 
measured by a condensation particle counter, (d) Classification of air mass based on 
back- trajectory analyses. Note in (a), the maximum concentration in the color scale is set 
at 20,000 cm” in order to show details of the size distribution evolution. 
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Figure 3. Eight clusters of 219 back- trajectories and their grouping into four back- 
trajectory classes. The back- trajectories were calculated for air parcels arriving at the 
SGP site 150 m AGL every three hours during the IOP. 
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Figure 4. Sub-micrometer aerosol volume concentration plotted again simultaneous 
measurements of CO mixing ratio for each air mass class. 
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Figure 5. Measurements during a nucleation event from 07:00 to 19:00 on May 7. (a) 
Dry aerosol number size distributions with the color indicating the number concentration, 
(b) SO? mixing ratio and number concentration of particles having diameters greater than 
10 nm. 
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Figure 6. Maximum SO 2 mixing ratio and particle number concentrations measured 
during the IOP as functions of wind direction. 
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Figure 7. Vertical profiles of measurements taken during an ascent from 15:40 to 16:12 
(local time) on May 25, 2003. (a) Relative humidity and potential temperature, (b) 
Aerosol scattering and absorption coefficients at 530 nm, adjusted to STP (0°C and 1 
atm), (c) Number concentration of particles having diameters greater than 10 nm, 
measured by a Condensation Particle Counter, concentration adjusted to STP, and (d) 
Aerosol number size distribution adjusted to STP. 
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Figure 8. Result of back- trajectory analysis for elevated aerosol layer observed above the 
SGP on May 25, 2003. 
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Figure 9. Aerosol number size distributions averaged over elevated layers during three 
flights and during the observation of plumes from Central America at the GIF. 
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Figure 10. Measurements from May 16, 12:00 to May 17, 15:00. (a) Ambient 
temperature and RH. (b) Sub-micrometer aerosol volume concentrations calculated from 
aerosol number size distributions, (c) Aerosol number size distributions, with the color 
indicating the number concentration, (d) Total particle number concentration ( D p > 1 0 nm) 
and precipitation. 
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Figure 11. (a) Aerosol size distributions at the beginning and the end of the first case for 
growth law analysis, and (b) derived particle growth rate and the least square fit to 
growth rates of condensational processes (Eqn. 5) and volume based reactions (Eqn. 6). 
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Image 1 . Extensive forest fires in Russia near the border with China shown in a picture 
taken by Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra 
Satellite on May 14, 2003 (From online MODIS rapid response system gallery, NASA). 
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